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ABSTRACT 

We report here a factor of 5.7 higher total CO flux in Arp 244 (the "Antennae" galaxies) than 
that previously accepted in the literature (thus a total molecular gas mass of 1.5xl0 10 M Q ), 
based on our fully sampled CO(l-O) observations at the NRAO 12m telescope. Currently, much 
of the understanding and modeling of the star formation in Arp 244 has been derived from a 
much lower molecular gas mass. It is imperative to reconsider all as the high molecular gas mass 
might provide sufficient fuel for ultraluminous extreme starburst in Arp 244 once the merging 
advances to late stage. 

Our observations show that the molecular gas peaks predominately in the disk-disk overlap 
region between the nuclei, similar to the far-infrared (FIR) and mid-infrared (MIR) emission. 
The bulk of the molecular gas is forming into stars with a normal star formation efficiency (SFE) 
Lir/M(H2) re 4.2Lo/M Q , same as that of giant molecular clouds in the Galactic disk. Additional 
supportive evidence is the extremely low fraction of the dense molecular gas in Arp 244, revealed 
by our detections of the HCN(l-O) emission, which traces the active star-forming gas at density 
>10 4 cm" 3 . 

Using the high-resolution BIMA plus the NRAO 12m telescope, full-synthesis CO(l-O) images, 
and our VLA continuum maps at 20cm, we estimate the local SFE indicated by the ratio map of 
the radio continuum to CO(l-O) emission, down to kpc scale. Remarkably, the local SFE stays 
roughly same over the bulk of the molecular gas distribution. Only some localized regions show 
the highest radio-to-CO ratios that we have identified as the sites of the most intense starbursts 
with SFE^20 L Q /M@ . Here we have assumed that the 20cm emission is a fairly good indicator 
of star formation down to kpc scale in Arp 244 because of the well-known correlation between the 
FIR and the radio continuum emission. These starburst regions are confined exclusively in the 
dusty patches seen in the HST optical images near the CO and FIR peaks where presumably the 
violent starbursts are heavily obscured. Nevertheless, recent large-scale star formation is going on 
throughout the system (e.g., concentrations of numerous super-star clusters and MIR "hotspots" ), 
yet the measured level is more suggestive of a moderate starburst (SFE^IO L Q /M©) or a weak 
to normal star formation (SFE^ 4 Lq/Mq), not necessarily occurring at the high concentrations 
of the molecular gas reservoir. The overall starburst from the bulk of the molecular gas is yet to 
be initiated as most of the gas further condenses into kpc scale in the final coalescence. 

Subject headings: galaxies: individual (NGC 4038/9) — galaxies: interactions — galaxies: starburst — 
infrared: galaxies — ISM: molecules — stars: formation 
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1. INTRODUCTION 

The "Antennae" galaxies (Arp 244, NGC 4038/39, 
VV 245), the nearest IR- luminous and perhaps the 
youngest prototypical galaxy-galaxy merger (the 
first in Toomre (1977) merger sequence), certainly 
reclaimed its fame from the recent releases of the 
HST/WFPC2 images (Whitmore et al. 1999) and 
the Chandra X-ray images (Fabbiano et al. 2000), 
and the observations at essentially all available 
wavelengths from radio to X-ray (e.g., Hummel & 
van der Hulst 1986; Neff, & Ulvestad 2000; Hib- 
bard, van der Hulst, & Barnes 2001 in preparation; 
Bushouse, Telesco, & Werner 1998; Nikola et al. 
1998; Mirabel et al. 1998; Vigroux et al. 1996; 
Read, Ponman, & Wolstencroft 1995; Fabbiano, 
Schweizer, & Mackie 1997). The HST images re- 
veal over one thousand bright young star clusters 
that are thought to have formed in recent bursts of 
star formation. Ha imaging also shows the most 
recent locations of the star-forming giant HII re- 
gions and their velocity fields as well (Rubin, Ford, 
& D'Odorico 1970; Amram et al. 1992; Whitmore 
et al. 1999). Soft X-ray (Fabbiano et al. 1997) 
and radio continuum (Hummel & van der Hulst 
1986) images may hint of previous star formation 
sites currently seen as supernova remnants. In ad- 
dition, mid-infrared (MIR) (Mirabel et al. 1998; 
Vigroux et al. 1996), far-infrared (FIR) (Evans, 
Harper, & Hclou 1997; R. Evans 1998, private 
communications; Bushouse et al. 1998) and sub- 
millimeter (Haas et al. 2000) images indicate that 
the most intense starburst takes place currently in 
an off-nucleus region that is inconspicuous at op- 
tical wavelengths. The star-forming molecular gas 
in Arp 244 is perhaps least understood, however; 
we show in this paper that the total molecular 
gas mass accepted in the literature for the last 
decade (Sanders & Mirabel 1985; Stanford et al. 
1990) has been underestimated by nearly a factor 
of 6. 

The importance of gas during galaxy-galaxy 
merging far exceeds its mass proportion, as 
demonstrated by the sophisticated numerical sim- 
ulations (e.g., Barnes & Hcrnquist 1996, 1998). 
At a distance of 20 Mpc (H = 75km s" 1 Mpc" 1 , 
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e.g., van der Hulst 1979; Mihos, Bothun, & 
Richstone 1993; Mirabel et al. 1998), the total 
IR luminosity (8-1000 ^m, Sanders & Mirabel 
1996) of Arp 244, measured from the IRAS four- 
band fluxes given in Soifer et al. (1989), is 
Lm = 6.2 x 10 10 L©. Using the most recent remea- 
sured total IRAS fluxes (as opposed to simply the 
"point source" in Soifer et al. 1989) in the Revised 
Bright Galaxy Sample (D.B. Sanders 2000, pri- 
vate communication; Sanders, Mazzarella, Kim, 
& Surace 2001 in preparation), the same IR lu- 
minosity is obtained. Therefore, strictly speak- 
ing Arp 244 is not a luminous infrared galaxy 
(LIG, Lir^10 11 L q ) unless a Virgocentric flow 
distance of 29.5 Mpc is used (Sanders et al. 
2000, but we use 20 Mpc throughout this pa- 
per), which leads to Lfir = 1-0 x 10 11 L Q and 
iiR = 1.3 x 10 11 L Q . A much smaller molecu- 
lar gas mass of 2.6 x 10 9 M Q (Sanders & Mirabel 
1985), currently accepted in the literature, leads to 
a higher global ratio of L m /M(H 2 ) ~ 24 L©/M Q . 
This ratio is often referred as the star formation 
efficiency (SFE) since the FIR emission, tracer of 
current star formation rate, has been normalized 
by the molecular gas mass available to make stars. 
With this high yield of young stars per unit molec- 
ular gas mass and a lower molecular gas content, 
most of the molecular gas will be depleted in ~ 10 8 
years, as noted by Sanders & Mirabel (1985) and 
Stanford et al. (1990). There would be no chance 
for Arp 244 to become an ultraluminous infrared 
galaxy (ULIG, Lir^10 12 L q ) in the late stage of 
the merging. Numerical modeling (Mihos et al. 
1993), based on this lower gas mass, indeed pre- 
dicted that Arp 244 would not join the rank of 
ULIGs. 

ULIGs, with the IR luminosity comparable to 
the bolometric luminosity of QSOs, are the most 
luminous galaxies in the local universe. They 
are believed to be powered mainly by starbursts 
(e.g., Smith, Lonsdale, & Lonsdale 1998; Gen- 
zel et al. 1998) taking place predominantly in 
the extreme starburst regions of a characteris- 
tic size of ~ 100 pc and Lir ~ 3 x 10 11 L Q 
(Downes & Solomon 1998), produced as a result 
of the merging of molecular gas-rich spiral galax- 
ies (Sanders & Mirabel 1996). Alternatively, the 
dust-enshrouded active galactic nuclei (AGN) may 
still be responsible for significant energy output 
in some ULIGs (Sanders et al. 1988; Veilleux, 
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Sanders, & Kim 1999; Sanders 1999). Now, for a 
total molecular gas mass of 1.5 x 1O 1O M , based on 
our fully sampled CO(l-O) observations obtained 
at the NRAO 6 12m, many previous conclusions 
and speculations about Arp 244 need to be revised. 
After all, the initial gas content, particularly the 
molecular gas — the fuel for star formation, will 
probably be the most important factor in deter- 
mining whether a merging pair of spiral galaxies 
reaches the peak of ultraluminous extreme star- 
burst phase since all ULIGs are still gas-rich with 
~ 10 10 M Q molecular gas mass (Solomon et al. 
1997), most likely concentrated in the kpc scale 
disks/rings around the merging nuclei (Downes 
& Solomon 1998; Scovillc, Yun, & Bryant 1997; 
Sakamoto et al. 1999; cf. Evans, Surace, & Maz- 
zarella 2000). 

With a molecular gas mass comparable to 
that of ULIGs and an early/intermediate merg- 
ing stage, Arp 244 is perhaps an example of what 
an ULIG might have looked like a few hundred 
million years ago. Arp 244 may be a snapshot 
in the evolution of a typical gas-rich merger into 
an ULIG system. Thus, it is crucial to under- 
stand how, where and when the starbursts initi- 
ate or have occurred in such an on-going merger. 
This becomes especially imperative in Arp 244 
given the advantage of its close-up distance and 
therefore the better linear resolution available. 
Nonetheless, the strongest starburst site revealed 
from the Infrared Space Observatory (ISO) MIR 
images (Vigroux et al. 1996; Mirabel et al. 1998; 
Xu et al. 2000) seems to be offset from the peak 
emission in all the 60, 100, and 160 /im FIR maps 
of the limited resolution, obtained with the Kuiper 
Airborne Observatory (KAO, Evans, Harper, & 
Helou 1997; Bushouse et al. 1998). In order to 
best locate the sites of the intense star formation, 
high-resolution imaging in the FIR is ultimately 
required. 

Since there is an excellent correlation between 
the FIR and the radio continuum emission (e.g., 
Helou, Soifer, & Rowan-Robinson 1985; Condon 
et al. 1990; Condon 1992; Xu et al. 1994; Marsh 
& Helou 1995), the FIR emission can thus be ap- 
proximately scaled according to the radio contin- 
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uum emission with the high resolution achievable 
by the Very Large Array (VLA) observations. We 
have therefore, obtained the VLA radio contin- 
uum images to compare with our high-resolution 
full-synthesis Berkeley-Illinois-Maryland Associa- 
tion (BIMA) plus the NRAO 12m CO maps. This 
is because the star formation occurs within gi- 
ant molecular clouds (GMCs) , especially the dense 
cores and starburst can be better characterized by 
the elevated SFE, i.e., the FIR to CO ratio, ap- 
proximated here by the radio-to-CO ratio. There- 
fore, local SFE across over the merging system can 
be approximately measured using the ratio map 
of the radio continuum to the CO emission to lo- 
cate the most intense starburst sites of the highest 
radio-to-CO ratios. 

We here briefly report our various observations 
in §2 detailing the single-dish CO mapping in 
Arp 244. The BIMA intcrfcromctry CO and the 
VLA 20cm continuum observations will be pre- 
sented elsewhere. The observational results and 
analysis are summarized and compared in §3. Sec- 
tion 4 discusses the starburst properties of Arp 244 
and the implications of our observations. Finally, 
we conclude our main results. 

2. OBSERVATIONS 

We obtained the NRAO 12m single-dish CO(l- 
0) map at Nyquist sampling (half-beam spacing 
~ 27. "5) in order to determine the true total 
CO extent and distribution and help add in the 
zero-spacing flux missed from the multifield BIMA 
CO(1-0) data cube. A total of more than 50 po- 
sitions have been observed in 1998 March and 20 
more positions in the outer edge of the map were 
further obtained in 1998 April to fully cover the 
entire CO emission region. We have essentially in- 
tegrated for at least 1.5 hr at every position in or- 
der to make sure that not only no further ^3c level 
CO emission is detected in the outmost edge of 
the map but all spectra obtained at different posi- 
tions have roughly comparable RMS noise level as 
well. A total of 73 positions have been mapped in 
the merging disks accumulating more than 110 hr 
useful integration time. Additional observations 
in 1999 June and November have been conducted 
remotely in the southeastern (SE) extension of the 
disk overlap region (where the much longer south- 
ern tidal tail starts) and at the tip of the southern 
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tail (Figure 1). This is to further integrate down 
the noise level to clearly show that CO is detected 
significantly far away from the merging disks, and 
to possibly search for CO emission at the loca- 
tion of the tidal dwarf galaxy several tens of kpc 
away from the merging disks. The typical inte- 
gration time at each of these selected positions is 
more than 3 hr. We have tried several positions 
at the tip of the southern tidal tail and deep in- 
tegrations in the possibly detected positions were 
further conducted in 2000 March and April, to- 
taling 27 hr usable integration time in this tidal 
feature area. 

We used the dual SIS 3mmhi receivers con- 
nected with both the two 256 x 2 MHz filter banks 
and the two 600 MHz spectrometers providing 
a velocity resolution of 2MHz <~ 5 kms -1 and 
a total velocity coverage of 1330 kms -1 . The 
system temperature T sys (SSB) was typically less 
than 400 K (on a scale) and the weather con- 
ditions were excellent throughout almost all ob- 
serving runs. Occasionally, exceptional weather 
with T sys <~ 200 K (for a low-decl. southern sky 
source at 115 GHz!) was seen in a few days. All 
observations were performed using a subreflector 
nutating at a chop rate of 1.25 Hz with a beam 
throw of ±3' plus a position switch (the so-called 
"BSP" mode) to achieve the flat baselines. This 
ensures that the off-source reference sky position 
is 6' away from the observed on-source pointing 
so that the telescope beam (FWHM~ 55") at the 
reference position is well outside the CO extent 
of Arp 244. Pointing and focus have been moni- 
tored frequently every 1-2 hr by observing nearby 
quasars 3C279 or 3C273 (and occasionally planets 
at the beginning of the observations). Uncertain- 
ties in positioning are typically ~ 5". Calibration 
with the standard chopper-wheel method was per- 
formed once after every the other 6 minute integra- 
tion scan and yielded an antenna temperature on a 
Tft scale. Further absolute flux calibrations in the 
antenna temperature scale have been performed 
by the repeated observations of the northern nu- 
cleus (our map center) during different observing 
sessions and by the observations of some nearest 
well-observed starburst galaxies. Comparing the 
observed line strengths and profiles, we found that 
the consistency is satisfactory and the difference is 
less than ~ 20%. 

HCN(l-O) was observed in 1997 April at only 



two locations: the northern nucleus (NGC 4038) 
and the CO peaks in the overlap region (Figure 1). 
This was part of another project to search for 
the HCN emission in the LIG mergers. Addi- 
tional integrations at these two locations were ob- 
tained in 1999 November as a further consistency 
check. The 12m telescope's FWHM beam at 89 
GHz (~ 72"), pointed at the nucleus of NGC 4038, 
presumably covers all HCN emission in NGC 4038, 
yet excludes emission from that of NGC 4039 and 
most of the disk overlap region. The second HCN 
beam covers not only all the CO peaks and the 
extended CO in the overlap but also most CO 
emission in NGC 4039. Therefore, the total HCN 
emission from Arp 244 can be roughly sampled 
by simply summing up these two beam measure- 
ments. T sys on a Tft scale is about 230 K with 
the dual SIS 3mmlo receivers and the same back 
ends as used in the CO observations. We have ac- 
cumulated a total integration time of nearly 5 and 
7 hr for NGC 4038 and the ovcrlap/NGC 4039, 
respectively. 

Our data reduction was performed using the 
CLASS/GILDAS package. Each individual scan 
was checked for spikes/dips or bad channels and 
curvatures in the baseline. Scans with the struc- 
tured baselines have been abandoned, bad chan- 
nels have been "repaired" by interpolation using 
the adjacent channels, and a linear baseline was 
then subtracted after removal of the spikes/dips 
for each accepted scan. All scans at each same 
position have then been summed to obtain an av- 
erage spectrum. The data cubes of the different 
velocity spacings were also created from the high- 
est velocity resolution (2MHz ~ 5 kms -1 ) and 
the smoothed (e.g., 4MHz ~ 10 kms -1 ) spectra 
maps. All spectra presented in this paper have 
been smoothed to 8 MHz (~ 21 kms -1 for CO 
and - 27 kms -1 for HCN). 

3. RESULTS AND ANALYSIS 
3.1. Previous CO(l-O) Observations 

The currently accepted total molecular gas 
mass in Arp 244 is based on a single point- 
ing observation with the 12m in the overlap re- 
gion. Sanders & Mirabel (1985) listed an in- 
tegrated intensity of 15.9 K kms -1 , or a total 
CO flux of only 556.5 Jy kms -1 (~ 35 Jy/K on 
a T£ scale), leading to a molecular gas mass of 
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2.6 xlO 9 Mq. The Owens Valley Radio Observa- 
tory (OVRO) millimeter array's two overlapping 
field (FWHM~ 65") CO imaging only recovered 
- 70% CO flux (Stanford et al. 1990) of the one- 
beam (FWHM~ 55") 12m measurement! The 
12m beam observed in Sanders & Mirabel (1985) 
only covered part of the overlap region, neither 
the nuclear regions of NGC 4038 and NGC 4039 
nor the extended structures throughout the entire 
merging disks (see Figures 1 and 2). Apparently, 
much of the CO is distributed over a much larger 
area than the 55" beam (Figure 2). Furthermore, 
the integrated CO line intensity in the old 12m 
observation was also underestimated owing to the 
insufficient velocity coverage (Figure 3, more in 
§3.2). 

Young et al. (1995) observed 6 positions in 
NGC 4038 and one position in NGC 4039 with 
the Five College Radio Astronomy Observatory 
(FCRAO) 14m (FWHM - 50") and derived a to- 
tal CO flux of ~ 2070 Jy km s -1 assuming an uni- 
form CO disk of radius - 0/6 for NGC 4038 and 
an exponential CO disk of a scale-length ~ 0/32 
for NGC 4039. Aalto et al. (1995) observed 3 
positions in Arp 244 with the Swedish European 
Submillimeter Telescope (SEST) 15m (FWHM ~ 
45"), i.e., the overlap region and the nuclear re- 
gions of NGC 4038 and NGC 4039. A sum of the 
integrated CO line intensities from the SEST ob- 
servations leads to a total CO flux close to what 
has been estimated by Young et al. (1995). Never- 
theless, the new OVRO three-field imaging (Wil- 
son et al. 2000) appears to only recover a total 
CO flux of 910 Jy kms" 1 , but this is already more 
than double that of the old OVRO map of Stan- 
ford et al. (1990). All these recent measurements 
of still very limited spatial coverage have already 
suggested that the total CO flux in Arp 244 can 
be a factor of a few larger than what has been 
previously accepted. 

3.2. New CO (1-0) Observations and the 
Total Molecular Gas Mass 

Figure 2 shows all the averaged CO spectra 
made with the 12m at a half-beam spacing, a to- 
tal of 73 spectra with more than 50 firm detec- 
tions. Only CO spectra at the outermost edge 
generally show non-detection or tentative detec- 
tions (noted in Table 1). Many CO spectra in 
the inner part of the map (the disk overlap and 



the two nuclear regions) have the integrated line 
intensities more than ~ 50% larger than that of 
the old 12m measurement of Sanders & Mirabel 
(1985). Observations in 1999 June at the same 
position as that of Sanders & Mirabel (1985) gave 
an integrated line intensity exactly twice larger, 
consistent with those mapped in the nearby posi- 
tions in 1998 March/April. The northern nucleus 
was always observed in all observing sessions giv- 
ing consistent results of about same integrated CO 
line intensity. Figure 3 presents the three spec- 
tra obtained near the CO peaks in the overlap 
region and compares with the original spectrum 
of Sanders & Mirabel (1985). Obviously, the old 
12m observation may have suffered severely from 
the poorly determined baseline owing to the lim- 
ited bandwidth, as the broadest CO spectra in the 
overlap region have FWZI ~ 600 kms -1 . But, this 
is probably insufficient to explain a factor of 2 dif- 
ference in the integrated CO line intensity. Ad- 
ditional errors in the telescope pointing and cali- 
bration together with the limited velocity coverage 
may have all contributed to the discrepancy. 

The total CO flux, from our fully-sampled 
CO observations (summarized in Table 1), is 
S CO dV = 3,172 ± 192 Jy kms -1 (7 C = 
90.6 ± 5.5 K kms -1 , using 35 Jy/K conversion 
for the T£ antenna temperature scale), a factor 
of 5.7 larger than what was initially reported by 
Sanders & Mirabel (1985). The total CO luminos- 
ity is thus Lqo = 0.3 x 10 10 K kms -1 pc 2 (T mb ), 
or a molecular gas mass (1.5±0.1) x 10 10 M (us- 
ing the standard CO-to-H 2 conversion factor of 
3.0xl0 20 H 2 cm- 2 (Kkms -1 ) -1 , or M(H 2 ) = 
4.78 x (L co /K kms -1 pc 2 ) Mq, applicable to 
GMCs in the Milky Way disk, e.g., Solomon & 
Barrett 1991; Young & Scoville 1991). The error 
estimate is purely from the measurement uncer- 
tainties (statistical errors from all measurements 
and the uncertainties in the CO detections in the 
outer edge), additional uncertainties up to ~ 20% 
should be expected from the calibration and point- 
ing etc. 

It is interesting to note that the new estimate 
of the dust mass in Arp 244, based upon the 
SCUBA 450 and 850 /xm measurements as well as 
the ISOPHOT observations at wavelength longer 
than 100 /xm, reveals a substantial amount of cold 
and warm dust, and leads to a total dust mass of 
M dust ~ 10 8 M Q (Haas et al. 2000). Thus, for 
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a total gas mass of M gas ~ M(H 2 ) + M(HI) ~ 
1.9 x 1O 1O M [M(HI) - 4 x 10 9 M Q , van der Hulst 
1979; Hibbard et al. 2001 in preparation], a gas- 
to-dust ratio of M gas /Mdust ~ 190 is obtained. 
This is basically the same gas-to-dust ratio estab- 
lished in our Galaxy and in some nearby spiral 
galaxies like NGC 891 (Alton ct al. 2000). This 
simple consistency may also imply that the stan- 
dard CO-to-H 2 conversion determined from GMCs 
in the Milky Way disk is likely applicable here and 
Arp 244 is indeed molecular gas-rich. Table 2 sum- 
marizes the global quantities of Arp 244. 

3.3. Global Molecular Gas Distribution 
and Kinematics 

Although the most molecular gas is concen- 
trated in the disk overlap region, the extended dis- 
tribution of the molecular gas spreads essentially 
all over the merging disks, much further beyond 
the field of view of the WFPC2 (Figure 2). There 
are also concentrations of the molecular gas in the 
two nuclear regions, particularly in the northern 
one. This is better indicated in the quite coarse 
channel map of 100 kms -1 width (Figures 4a- 
c) . In general, there are good correspondences be- 
tween the CO distribution and the optical mor- 
phology, except for the overlap region where the 
CO emission appears to be much more prominent 
than the optical light. Figure 4f further compares 
the CO integrated intensity contours (using all 74 
position CO observations) with the optical disks. 
All these (Figures 4a-4f) clearly show that the 
molecular gas extends throughout the system and 
explain why the old 12m single beam measurement 
underestimated the total CO emission by a large 
factor. Obviously, much of the CO over a broad 
velocity range is simply distributed over several 
arcminutes, corresponds to a linear CO extent of 
possibly beyond ^20 kpc. 

The deep integration at the SE edge of the over- 
lap region, where the southern tidal tail begins, 
confirms the clear detections of the weak CO emis- 
sion from this region (Figure 5a). The central off- 
set of the spectra shown in Figure 5a is (110", 
—82.5") relative to the northern nucleus, two full 
beams further out from the peak emission in the 
overlap region. CO emission is also detected, for 
example, at offsets of (55", 55") and (81", 27.5") 
relative to the northern nucleus, and in a few lo- 
cations west of the two nuclear regions at one and 



half beams from the nuclei (Figure 2 and Table 1). 
These locations are so far away from the CO emis- 
sion peaks and their immediately associated CO 
extensions (also see the new OVRO map of Wil- 
son et al. (2000)), suggesting that there are weak 
extended CO features across over the entire op- 
tical disks. Given the relatively quite strong CO 
emission and the large distance from the CO emis- 
sion peaks, it is unlikely that both the pointing 
errors, as large as 10" in some rare occasion, and 
the telescope beam pattern (sidelobes) can pro- 
duce such relatively strong CO emission that we 
have detected near the map edge. This certainly 
points out the importance of the complete cov- 
erage of the merging systems in the single beam 
observations since they may systematically under- 
estimate the total CO content in the cases where 
the merging disks are larger than the single beam. 
A substantial amount of the molecular gas may be 
distributed at large galactic radii, especially in the 
young, early and intermediate stage mergers like 
the "Antennae" galaxies. 

The overall velocity spread of the molecular gas 
is more than 600 km s _1 across over the two merg- 
ing disks. The narrowest velocity ranges are ob- 
served in the nuclear region of the northern face-on 
galaxy NGC 4038 (FWZI - 300 kms" 1 , FWHM - 
100 kins" 1 ), whereas a broader velocity spread 
(FWZI > 500 kins" 1 , FWHM - 230 kms" 1 ) 
can be seen in the southern inclined galaxy 
NGC 4039. The broadest velocity components 
(FWZI - 600 kms" 1 , FWHM - 300 kms" 1 ) are 
particularly prominent south of NGC 4039 even 
though the CO emission is quite weak (Figure 2 
and Nos. 6, 9, 17, 27, & 29 in Table 1). It is 
likely that the interferometric observations have 
almost entirely resolved and thus missed these 
weak and broad velocity features. The velocity 
extent around the peak CO emission in the over- 
lap has a similarly broad velocity range. Also, 
the velocity spread of most spectra in the over- 
lap region lies between the two extremes observed 
in the two galaxies. This is obvious since some 
spectra contain the CO contribution from either 
NGC 4038 or NGC 4039 or both. But, clearly 
most spectra at a full beam (or more) eastern of 
the nuclei are basically free from any contribution 
of the nuclear CO emission. We also notice that, 
although the difference is small (~ 50 kms -1 ), 
the systemic velocity in the overlap region is lower 
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than that of the entire system whereas the sys- 
temic velocity south of NGC 4039 is higher. This 
perhaps reminisces the disk rotation of the molec- 
ular gas in the more inclined galaxy NGC 4039. 
Nonetheless, the average velocities of the molecu- 
lar gas in the nuclear regions of both galaxies are 
about same as the systemic velocity of the entire 
system. 

The detailed study of the molecular gas distri- 
bution and kinematics is difficult given the lim- 
ited resolution of the 12m telescope beam, which 
will be thoroughly described in the presentation 
of the BIMA + NRAO 12m full- synthesis obser- 
vations, providing a much improved spatial res- 
olution of nearly 10 times better. Although the 
enormous amount of molecular gas in the overlap 
region must originate from both spiral galaxy pro- 
genitors, the kinematics and the proximity of the 
molecular gas in the overlap, combined with the 
similar velocity spread of the molecular gas motion 
as that in the southern galaxy NGC 4039, seem to 
suggest that the most molecular gas in the overlap 
could be originated from the southern progenitor. 
The molecular gas content of the two progenitors 
might be initially comparable and both were gas- 
rich. However, the VLA HI observations appear 
to show that the progenitor of NGC 4039 was gas 
poor since the most HI in this system is currently 
in the southern tail extending from NGC 4038 (Hi- 
bbard et al. 2001 in preparation). 

3.3.1. CO in the Southern Tidal Dwarf Galaxy 

We have also searched for the faint possible CO 
emission at several positions in the southern tidal 
dwarf galaxy (Schweizer 1978; Mirabel et al. 1992) 
at the tip of the tidal tail. A possible weak CO de- 
tection from the beam at position RA=12:01:26.0, 
Dec=- 19:00:37.5 (J2000) appears to be at the 
~ 4cr level (an integration of more than 12 hr; 
Fig. 5b). An average of the summed spectra of all 
observations at the 5 different positions we have 
searched (accumulated an integration time of more 
than 27 hr), appears to confirm the weak CO de- 
tection at the ~ 5ct level. All these observed po- 
sitions are selected to be around the regions of 
the highest HI column density, where the VLA HI 
observations (Hibbard et al. 2000) have revealed 
some large concentrations of the atomic gas, dis- 
tributed extendedly around the region of the tidal 
dwarf galaxy and along the southern tidal tail. 



The weak CO detections (Figure 5b) have also 
shown the same velocity range as the atomic gas 
velocity spread, revealed from the VLA HI channel 
maps, further indicating that the weak CO emis- 
sion is likely real. Yet the much deeper integration 
and the more observing positions in this region are 
still required in order to firmly present the CO de- 
tections at the very high significant levels, and to 
possibly obtain a better understanding of the weak 
CO distribution and kinematics in a tidal dwarf 
galaxy. 

The integrated CO line intensity from the de- 
tected position is Iqo = 0.37 ± 0.09 K kms -1 . 
The average spectrum of the sum from all other 
4 locations (not shown here) appears to give a 
similar integrated CO line intensity Iqo — 0.35 ± 
0.12 K kms -1 , while the average spectrum of all 
observations gives an integrated CO line intensity 
-^CO — 0.35 ± 0.07 K kms -1 . Earlier observa- 
tions by Smith & Higdon (1994) failed to detect 
any CO emission, but our observations here are 
much more sensitive. Using the standard CO-to- 
H2 conversion factor, we therefore obtain a mini- 
mum molecular gas mass for the tidal dwarf galaxy 
~ 0.61 x 10 8 M© if only the detected position is 
considered. It is likely that the weak CO emis- 
sion exists in other locations in this region and 
the total molecular gas mass for the entire tidal 
dwarf galaxy can be a factor of several larger. In- 
deed, the summed average spectrum of all the 5 
position observations seems to suggest a molecu- 
lar gas mass ^2x 10 8 M Q . Anyway, this molecu- 
lar gas mass appears to be larger than that of the 
molecular complexes detected in the M81 group 
(Brouillet, Hcnkel, & Baudry 1992; Walter & Hei- 
thausen 1999), yet comparable to the molecular 
gas mass detected in the other two tidal dwarf 
galaxies (Braine et al. 2000). It is possible that 
most of the molecular gas might have come di- 
rectly from the conversion of HI into H2 . This is 
similar to the case in the two tidal dwarf galaxies 
where the molecular gas was detected in regions of 
the highest HI concentrations (Braine et al. 2000), 
and is also related to the case in the cold intra- 
group medium where huge HI concentrations and 
the CO detection are both associated with the in- 
tragroup starburst (Gao & Xu 2000) in the famous 
compact galaxy group Stephan's Quintet. 
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3.4. HCN(l-O) Detections and the Dense 
Molecular Gas 

The high dipole-moment molecules like HCN, 
which traces the dense molecular gas at density 
£ 10 4 cm -3 , are better tracers of the star-forming 
gas than that of CO (Solomon et al. 1992; Gao 
1996; Gao & Solomon 2000a, 2000b). We have 
detected the; weak HCN(l-O) emission from the 
two positions, one at the nucleus of NGC 4038 
and the other covering both the disk overlap and 
the nuclear region of NGC 4039 (FWHM~ 72" 
at 89 GHz, the two thick circles in Fig. 1). Fig- 
ure 5c shows the two HCN spectra and the total 
HCN luminosity can be estimated by the sum of 
the two beam measurements. This is because, in 
nearby galaxies where HCN maps exist, the dense 
molecular gas tends to be concentrated in the in- 
nermost disks of the highest density regions (Gao 
1996; Heifer & Blitz 1997a). Thus it is unlikely 
that any significant HCN emission in Arp 244 still 
exists outside the coverage of the two HCN beams. 
We estimated the total HCN flux to be SucNdV = 
39 ±6 Jy kms" 1 (7 H cn = 1.2 ±0.2 K kms" 1 , us- 
ing a ~ 32 Jy/K conversion) and the HCN lumi- 
nosity Lhcn = 0.7 x 10 8 K kms" 1 pc 2 . There- 
fore, the luminosity ratio of £hcn/£co = 0.02 
implies only £5% of the molecular gas is at high 
density of £10 4 cm" 3 in Arp 244 (Gao 1996; Gao 
& Solomon 2000a). This is surprisingly low since 
ULIGs can have as high as Lhcn/^co = 0.25 and 
up to about half of the molecular gas at density 
£l0 4 cm" 3 (Solomon ct al. 1992; Gao 1996; Gao 
& Solomon 2000a). The ratio of £hcn/£co in 
Arp 244 is actually just comparable to that of the 
Milky Way disk (Heifer & Blitz 1997b), even lower 
than that of some "normal" spiral galaxies (Gao & 
Solomon 2000a) . The fraction of the dense molec- 
ular gas in Arp 244 is apparently the lowest among 
all LIGs observed so far. 

Although most of the CO emission (£ 60 %) 
is from the overlap region, the HCN emission in 
the overlap appears to be less than half since only 
about half of the HCN emission is originated from 
both NGC 4039 and the overlap. The HCN de- 
tection from the overlap region and NGC 4039 
has only a comparable integrated line intensity, 
but weaker in the antenna temperature as com- 
pared with the detection in NGC 4038. The two 
HCN spectra differ clearly in the velocity spread: 
a narrow velocity space in the HCN emission of 



NGC 4038 whereas a much broader velocity spread 
in NGC 4039 and the overlap region. Thus, the 
HCN spectra basically indicate that the dense 
molecular gas could have the similar gas kinemat- 
ics as that of the total molecular gas revealed from 
the CO emission even though their spatial distri- 
butions are probably totally different. 

The low fraction of the dense molecular gas in 
Arp 244, especially in the overlap region, is con- 
sistent with the fact that the bulk of the molec- 
ular gas in the overlap is extensively distributed 
both spatially and kinematically (Figs. 2, 3 and 
4). Although there are several CO peaks in the 
high resolution synthesis maps, none of them has 
the extremely high concentrations of molecular gas 
as those found in ULIGs. The highest gas sur- 
face density in the northern nucleus is ~ 2.5 x 10 3 
M Q pc~ 2 (resolved on a linear scale of ~ 0.5 kpc 
as in both the BIMA plus NRAO 12m and the 
new OVRO CO maps), order of magnitude smaller 
than those in ULIGs, whereas the gas surface den- 
sity in the southern nuclear region is lower by 
at least a factor of 3 and comparable to that of 
the CO peaks in the overlap. Our low-resolution 
(55" ~ 5 kpc) 12m map reveals a CO peak in the 
disk overlap region since most of the molecular gas 
is concentrated there. The peak molecular gas col- 
umn density N(Ha) is close to 10 22 cm~ 2 , i.e., an 
average gas surface density (over ~ 5 kpc scale) of 
about 190 M Q pc -2 . This low gas surface density 
together with the low HCN-to-CO line luminos- 
ity ratio clearly marks the unique gas properties 
in Arp 244, in sharp contrast with those in most 
LIGs/ULIGs. 

3.5. Local Star Formation Efficiency and 
the Radio-to-CO Ratio Map 

There is apparently a correlation between the 
radio continuum and the CO emission as even so 
indicated in the old OVRO CO map (Hummel & 
van dcr Hulst 1986; Stanford et al. 1990). Partic- 
ularly, most of the radio continuum emission is ex- 
tensively concentrated in the overlap region with a 
very similar morphology to that of the CO. The di- 
rect comparison of the KAO 60 fim maps (Evans 
ct al. 1997) with our full-synthesis CO images 
clearly reveals that the average FIR to CO ratio 
in the overlap region is higher than the global ratio 
since most of the FIR emission (£75%) is from the 
overlap region, whereas only £60% CO emission 
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is from the same region. Therefore, the average 
SFE in the overlap can be a bit higher than the 
global value (which is Lir/M(H 2 ) = 4.2 L Q /M Q , 
or 5.3 Lq /Mq if a higher 100 fj,m KAO flux of 
Bushousc et al. (1998) is used), still comparable 
to that of GMCs in the Galactic disk and an order 
of magnitude lower than that of ULIGs (Sanders & 
Mirabel 1996; Solomon ct al. 1997). On the other 
hand, the average SFE in the nuclear regions is 
expected to be lower than the global value since 
the two nuclear regions have extremely weak FIR 
emission, whereas the CO is highly concentrated, 
particularly in NGC 4038. Moreover, the radio 
continuum emission has indeed an excellent cor- 
respondence with the FIR maps, just as expected 
from the well-known correlation between the radio 
continuum and the FIR emission (e.g., Helou et al. 
1985; Marsh & Helou 1995). 

In order to further quantify the local star for- 
mation properties and better identify the sites of 
the most intense starbursts with the highest SFE, 
we have produced the ratio map of the VLA 20cm 
radio continuum to the CO(l-O) emission. Since 
the tight FIR/radio correlation is also valid on kpc 
scale in galaxies (e.g., Marsh & Helou 1995; Lu 
ct al. 1996), the radio continuum emission can be 
used as a tracer of the recent star formation in 
galaxies as well (e.g., Condon et al. 1990, 1991; 
Condon 1992). We here use the radio continuum 
maps to roughly indicate the star formation sites. 
Thus, a high resolution indicator of the FIR emis- 
sion can be inferred by scaling the FIR according 
to the radio continuum emission. Comparing the 
radio continuum maps with the detailed molecular 
gas distribution from the CO imaging, we can ob- 
tain the radio-to-CO ratio map to reveal the local 
SFE, which is defined as the local ratio of the star 
formation rate to the molecular gas mass. Both 
the VLA 20cm continuum and the BIMA plus the 
NRAO 12m full-synthesis CO maps have about 
same resolution, thus a direct division can be per- 
formed to obtain the ratio map. 

The radio continuum to the CO ratios have 
been plotted as contours and compared with the 
HST/WFPC2 (Fig. 6) and the VLA radio contin- 
uum images (Fig. 7). Apparently, a rather smooth 
distribution of the radio-to-CO ratio is observed 
across over most of the merging disks. The sites 
of the highest ratios, ~ 2.5 and 5 times above 
the average, have been indicated in Figure 6 as 



black and white contours, respectively, and they 
coincide well with the dusty patches across over 
the HST/WFPC2 image. The sites of the high- 
est radio-to-CO ratios are exclusively localized in 
the most prominent dusty regions in the overlap, 
which are near the CO peaks and probably the 
FIR emission peaks. The correspondent highest 
SFE is therefore ^,20 L Q /M Q since the global av- 
erage SFE is 4.2 L Q /M Q . Interestingly, almost 
all LIGs/ULIGs have SFE > 20 L Q /M (Sanders 
& Mirabel 1996), although some early stage pre- 
merging LIGs with SFE ~ 10 L Q /M Q and lower 
do exist (e.g., Arp 302, Lo, Gao, & Gruendl 1997). 

Consistent with what can be roughly expected 
from a direct comparison of the FIR map at 60 ^m 
with the CO images in the nuclear regions, both 
nuclei have only a low radio-to-CO ratio, about 
the average or lower (Fig. 7). Thus the SFE in the 
nuclear regions is about 4 L©/M Q . The ISO MIR 
"hotspots" have higher ratios than the average, so 
do the extended overlap region and the western 
star- forming loop in NGC 4038, in addition to a 
few spots around the nuclear regions (especially 
the circumnuclear region of NGC 4038, Fig. 6). 
But these are a factor of 2 lower than the high- 
est ratios and are only of moderate starburst sites 
with a SFE ^10 L Q /M Q , which is actually the 
typical SFE value for the local starburst galaxies 
(Sanders & Mirabel 1996). 

4. DISCUSSION 

Both observational and theoretical studies 
have demonstrated that the ultraluminous ex- 
treme starburst phase is most likely achieved in 
the molecular gas-rich mergers (e.g., Sanders & 
Mirabel 1996; Mihos & Hernquist 1996; Downes 
& Solomon 1998; Gao & Solomon 1999) where 
starbursts may proceed through the formation 
of numerous super-star clusters (e.g., Whitmorc 
et al. 1999; Surace et al. 1998; Scoville et al. 
2000). Multiwavclcngth observations of Arp 244 
and comparison with the early stage mergers ob- 
served by Gao et al. (1999) are therefore espe- 
cially crucial to track how and when the dominant 
sources of star formation transform from within 
the disks of the two gas-rich spirals to the disk-disk 
overlap regions between the two galaxies. More- 
over, comparative study with CO observations of 
the late stage mergers (Downes & Solomon 1998; 
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Bryant & Scoville 1999) are also equally impor- 
tant for understanding how the starbursts con- 
tinue to evolve from taking place predominantly 
in the overlap regions to occurring mainly in nu- 
clear regions when the merging advances to ULIG 
phase, with extraordinary nuclear concentrations 
of dense gas. Intermediate mergers like Arp 244 
serve apparently as an important link in the merg- 
ing process between a pair of gas-rich premergers 
and the merged double-nucleus ULIG. We here 
discuss only the starburst properties of Arp 244 
related to our 12m observations. 

4.1. Ultraluminous Extreme Starburst 
Potential 

The CO content is decreasing as merging pro- 
gresses, indicating the depiction of molecular gas 
due to merger-induced starburst (Gao & Solomon 
1999). Using the scaling of M(H 2 ) - (Gao 
& Solomon 1999, valid for mergers at the early 
and intermediate stages, with the projected nu- 
clear separation 20 > 5 scp > 2 kpc), the total 
molecular gas will decrease roughly by more than 
a factor of 2 when Arp 244 reaches S sep ~2 kpc, 
prior to or about entering into the late merging 
stage. It appears that all ULIGs observed so far 
have ~ 10 10 Mq of molecular gas mass or more 
without exception, even after considering up to a 
factor of 5 reduced CO-to-H2 conversion (Solomon 
et al. 1997). The ultraluminous starburst phase 
of Arp 244 is thus possibly reachable in the ad- 
vanced stage since only ~ 10 10 M Q molecular 
gas will be consumed by the newborn stars in 
the next ~ 10 8 yr of the merging process. There 
will still be abundant molecular gas of close to 
~ 10 10 M© available, not to mention the addi- 
tional amount of atomic gas in the merging disks 
and in the long tidal tails (Hibbard et al. 2001 in 
preparation) that may eventually fall back to the 
merging disks (Hibbard et al. 1994). The atomic 
gas could be an additional gas reservoir for possi- 
ble conversion into the molecular phase, because 
there appears to be some evidence that ULIGs 
tend to have the highest ratio of M(H 2 )/M(HI) 
(e.g., Mirabel & Sanders 1988, 1989). 

In general, the gas surface density increases 
from orders of a few times 10 2 M Q pc _2 to a 
few times 10 3 M Q pc~ 2 , as merging progresses 
from early to intermediate stages (Gao et al. 
1999) (unresolved on a scale of 1-2 kpc), whereas 



advanced ULIG mergers such as Arp 220 have 
gas surface density typically greater than 10 4 
M pc~ 2 (Scoville et al. 1997; Downes & Solomon 
1998; Sakamoto et al. 1999). Although the total 
molecular gas content of Arp 244 is comparable 
to that of ULIGs, the gas surface density is still 
orders of magnitude lower than that of ULIGs, 
e.g., in the overlap, ~ 10 2 M Q pc~ 2 on a scale 
of 5 kpc and <~ 10 3 M Q pc~ 2 on a resolved scale 
of sub-kpc. Therefore, only a small fraction of 
the molecular gas in Arp 244 is actually at a high 
enough gas surface density for a high SFE, as 
formulated in the Schmidt law (e.g., Kennicutt 
1998). The overall gas density is low to maintain 
a low SFE across over the entire system, except 
for some localized starburst regions (Fig. 6), just 
as our HCN observations have revealed that little 
dense molecular gas available currently to power 
the extreme starbursts. 

Ultraluminous extreme starbursts require not 
only a large quantity of molecular gas but also 
a high gas concentration, particularly the nuclear 
gas concentration, so that the bulk of the molec- 
ular gas is at a high density. The large quantity 
of pre-existing molecular ISM is the sufficient fuel 
for Arp 244 to ultimately reach the onset of the 
ultraluminous starburst phase, as the merging pro- 
ceeds to the advanced stage, and the bulk of the 
molecular gas becomes highly condensed. Unlike 
most LIGs/ULIGs, which show the dominant nu- 
clear MIR and perhaps FIR emission (Hwang et al. 
1999; Soifer et al. 2000; Bushouse et al. 1998; Xu 
et al. 2000), Arp 244 has only moderate nuclear 
emission in MIR (Mirabel et al. 1998; Xu et al. 
2000), and the KAO 60 /im emission from the 
two nuclei is extremely weak (Evans et al. 1997). 
Although there is no clear indication of a large 
stellar bulge in either of the galaxies, the concen- 
trations of the molecular gas as well as the stars in 
the nuclear regions in Arp 244 perhaps also draw 
some analogies to the case of galaxy-galaxy merger 
model with stellar bulges in the progenitors (Mi- 
hos & Hernquist 1996). In reality, Arp 244 is prob- 
ably more in between the two extreme cases mod- 
eled by Mihos & Hernquist (1996). According to 
the models (see Figure 5 in Mihos & Hernquist 
1996), the first peak of starburst phase has re- 
cently passed in Arp 244, shortly after the pair's 
first closest approach, perhaps as evidenced by nu- 
merous young star clusters in the HST/WFPC2 
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images. Thus, Arp 244 could have just phased 
out the bona fide starburst LIG stage. As more 
and more gas continues to accumulate in both 
the overlap and the two nuclear regions, the fur- 
ther nuclear gas concentrations and the continuous 
star formation can build stronger bulges that can 
help stabilize the disks against the bar formation, 
lower the SFE levels during the merging, and leave 
enough gas for the strongest final starburst when 
the two galaxies eventually coalesce, with almost 
all the molecular gas collapsed into a kpc double- 
nucleus region. More than an order of magnitude 
increase in the SFE, combined with an abundant 
molecular gas supply of ~ 10 10 M Q , will be just 
sufficient to make Arp 244 into Lir ~ 10 12 L Q . 
Interestingly, there also appears to be some quite 
widely separated ULIGs (5 sop > 10 kpc), which 
might still be in their early stage of merging, 
but may have passed their first closest impacts, 
and may have just been experiencing their first 
strongest starbursts (Murphy 2000). 

In summary, a preexisting abundant molecular 
gas content is a necessary for there to be a possibil- 
ity of a merging spiral pair reaching the ultralumi- 
nous starburst phase. Arp 244, with a large quan- 
tity of molecular gas available at low density, a low 
SFE, a moderate FIR luminosity or star forma- 
tion rate, and yet a relatively early/intermediate 
stage of merging, has the potential of producing 
an ultraluminous extreme starburst in late stage 
of merging, when the molecular gas complexes in 
the overlap region eventually merge with the nu- 
clear gas concentrations and finally collapse into a 
kpc-scale double-nucleus gas disk. 

4.2. Starbursts and the Current Intense 
Star Formation Sites 

The well-known correlation between the FIR 
thermal dust emission and the radio continuum 
synchrotron emission, unexpected as they are ap- 
parently two unrelated physical mechanisms, is 
probably the tightest relation known among the 
global quantities of galaxies (e.g., Helou ct al. 
1985; Condon 1992; Xu ct al. 1994). The mor- 
phologies of the two emission are also similar 
and the correlation between them appears to hold 
down to kpc scale within individual galaxies (e.g., 
Bicay, Condon, & Helou 1989; Marsh & Helou 
1995; Lu et al. 1996). Based on these, we approx- 
imately scaled the local FIR emission according to 



the radio continuum emission, and compared with 
the molecular gas distribution to measure the lo- 
cal SFE. We have clearly identified that the most 
intense bursts of star formation (the sites of the 
highest SFE, i.e., the highest radio-to-CO ratios, 
Figure 6), currently being observed, are not in the 
concentrations of the optically revealed star clus- 
ters, or even the peaks of the MIR emission, nor 
necessarily at the exact locations of the CO and 
20cm radio continuum emission peaks. Instead, 
the strongest starbursts appear to be confined in 
small regions near the CO, radio continuum, and 
FIR emission peaks in the overlap region. The 
bulk of the molecular gas across over the entire 
system is making stars with a rather modest or 
"normal" SFE (Fig. 7). 

Let us emphasize what we have learned from 
our observations and the various other observa- 
tions. Starbursts have apparently been going on or 
happened in some regions in Arp 244 as evidenced 
by super-star clusters (Whitmore et al. 1999) and 
the MIR "hotspots" in the overlap (Mirabel et al. 
1998). Yet, most of the FIR emission (as well as 
the MIR), which dominates the total energy out- 
put, comes from the overlap region (Evans et al. 
1997; Bushouse et al. 1998), where most molec- 
ular gas resides, rather than from the super-star 
cluster concentrations. Although Arp 244 has long 
been claimed as the nearest archetypal starburst 
merger, our observations show that the molecu- 
lar gas is rather extensively distributed over the 
entire merging disks, though most is still in the 
overlap, and the bulk of the molecular gas is only 
making stars at a low SFE. The entire system is 
presently not undergoing a global burst of star for- 
mation, even though some localized starbursts are 
occurring in the overlap, and a global starburst 
could have just peaked recently, right after the first 
closest impact responsible for the formation of the 
tidal tails. 

Whitmore et al. (1999) found that the star clus- 
ters at the edge of the dusty overlap region appear 
to be the youngest, with ages ,$5 Myr, while the 
star clusters in the western loop in NGC 4038 ap- 
pear to be 5-10 Myr old. Indeed, these are the 
sites of the most intense and moderate starbursts 
respectively as revealed by the radio-to-CO ratio 
map in Figure 6. There might be many youngest 
star clusters totally obscured by the dust at the 
most intense starburst sites, contributing signifi- 
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cantly to the dominant energy output in the over- 
lap region. On the other hand, many star clusters 
in the northeastern that appear to be ~ 100 Myr 
old, and even older star clusters across over the 
system, are not the current starburst sites, con- 
tributing little to the total FIR emission. These 
old star clusters are probably related to the recent 
large-scale active star formation happened after 
the first closest encounter. 

The excellent agreement among the sites of the 
highest radio-to-CO ratios, the dark patches in 
the HST/WFPC2 images, and the youngest ages 
of the star clusters, is striking. Apparently both 
the most intense localized starbursts and the over- 
all large-scale star formation are occurring in the 
overlap region, where all the CO, [C u] line (Nikola 
ct al. 1998), FIR, and radio continuum show 
high concentrations. But, the confined intense 
starburst sites, which are heavily obscured in the 
dust, have only SFE ^,20 L©/M Q , just reaching 
the typical SFE level for LIGs/ULIGs. On the 
other hand, the average SFE in the overlap is just 
slightly larger than 4 L Q /M , still comparable 
with that of the Galactic disk GMCs. 

One caveat is that Arp 244 may be twice lu- 
minous in 100 /im, as measured by the KAO 
(compared with that of the IRAS measurement, 
Bushouse et al. 1998), although the 60 /im mea- 
surements agree (Evans et al. 1997). This dis- 
crepancy needs future FIR observations to re- 
solve. More importantly, it is perfectly possible 
that the standard CO-to-H2 conversion factor can 
be an overestimate for the total molecular gas 
mass. Therefore, it is likely that the global SFE 
in Arp 244 can be up to more than a factor of 2 
larger than that of the GMCs in the Milky Way 
disk. This may indeed suggest that an enhanced 
global SFE proceeds throughout the entire merg- 
ing disks while some confined strongest starbursts 
have a highest SFE much larger than 20 L Q /M Q , 
typical for ULIGs. In any case, the highest radio- 
to-CO ratio sites could actually be bona fide sites 
of the current starbursts with an elevated SFE, 
and may contribute a significant fraction to the 
total FIR emission, even though most of the FIR 
emission may still come from the regions of a nor- 
mal or modest SFE. 

Using the [C n] /CO (1-0) line ratio to distin- 
guish between starburst activity in galaxies and 
more quiescent regions (Stacey et al. 1991), Nikola 



ct al. (1998) concluded that there is no strong 
starburst activity taking place in Arp 244 on a 
scale of ~ 5 kpc, the resolution of their data. 
This is roughly in agreement with our results since 
starburst sites are only identified on a scale of 
<1 kpc (Fig. 6). Nikola et al. (1998) further ar- 
gued that most of the [C n] emission might have 
come from the confined active star forming regions 
surrounded with the more quiescent GMCs. Our 
results appear to support this claim if indeed the 
sites of the intense starbursts, which are identi- 
fied from the high radio-to-CO ratios (Fig. 6), are 
producing the bulk of the [C n] emission. Ad- 
ditional supports for the confined starbursts in 
small regions in the overlap are the recent ISO- 
SWS measurements (Kunze et al. 1996) and the 
ISO-LWS observations, including the detection of 
Br7 knots in the overlap interaction zone (Fischer 
et al. 1996). 

The rather extended moderate starbursts also 
coexist closely with the confined, most intense 
starbursts in the overlap region. Other interest- 
ing sites of moderate starbursts are the western 
loop/ring structure, and part of the £ kpc scale 
circumnuclcar regions (mainly in NGC 4038), but 
not the nuclei themselves (Figs. 6 & 7), with nearly 
3 times higher SFE than the average. The west- 
ern loop coincides with the molecular ring we have 
mapped from the BIMA plus the NRAO 12m full- 
synthesis CO image, which was completely ab- 
sent from the old OVRO CO map (Stanford et al. 
1990), and was only partially revealed in the new 
OVRO CO map (Wilson ct al. 2000). It appears 
that the starbursts occur at a moderate level in the 
southeast side of the overlap region, proceed pro- 
gressively with increased intensity towards north- 
west across over the overlapping molecular con- 
centrations, and produce the most vigorous star- 
burst in the northern and western edges of the 
huge molecular gas agglomerations in the overlap 
interaction zone (Fig. 7). 

Moderate [C n] emission is also detected in the 
western loop (Nikola et al. 1998), which also in- 
dicates the moderate star-forming activity there. 
This is again consistent with our results obtained 
from the radio-to-CO ratios. Moreover, the [C n] 
emission peak appears to be slightly offset from 
our single-dish CO peak (Figures 2 & 4f), which 
is resolved into several CO peaks by the interfer- 
ometers. The strongest [C n] emission seems to 
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be roughly coincident with the sites of the high- 
est SFE, just north of the CO peak in the NRAO 
12m map. Our single-dish CO map also appears 
to be roughly peaked at same position as that of 
the KAO 60 (Evans et al. 1997), 100 and 160 {im 
maps (Bushouse et al. 1998). The strongest MIR 
peak in the overlap region that is inconspicuous at 
optical wavelength (Mirabel et al. 1998; Xu et al. 
2000), however, may correspond only to the south- 
ernmost CO peak in the overlap (Wilson et al. 
2000), which is not one of the strongest starburst 
sites and has only a moderate SFE ^,10 L Q /M 
(Figures 6 & 7). Therefore, the strongest peaks of 
the CO, radio continuum, and FIR, as well as the 
MIR peaks, arc not necessarily the exact sites of 
the highest SFE, or the most intense starbursts. 
The highest radio-to-CO ratio peaks, which may 
reveal the sites of the most vigorous star for- 
mation with the highest SFE, and the strongest 
[C n] emission peaks, which arise mainly from 
the photodissociation regions, appear to be the 
best probes of the most intense starbursts, which 
might be totally obscured in the optical/near-IR 
and even in the MIR regime. 

4.3. Star Formation from the Multiwave- 
length Observations 

In light of the numerical models for mergers of 
Arp 244 alike (Toomre & Toomrc 1972; Barnes 
1988; Mihos et al. 1993; Mihos & Hernquist 1996; 
Barnes & Hernquist 1996, 1998) and the multi- 
wavelength observations of Arp 244, comparison 
of the various observations and statistical studies 
of LIG mergers with those models can be fruitful. 
We here try to discuss qualitatively the onset of 
the starbursts, the fate of the starbursts, and the 
possible ultraluminous extreme starburst phase in 
Arp 244, utilizing these observational and theoret- 
ical results. 

ROSAT high-resolution imager (HRI) X-ray 
map (Fabbiano et al. 1997) and the high- 
resolution and high-sensitivity Chandra X-ray im- 
ages from the Chandra news release (Fabbiano 
et al. 2000), after the submission of this paper, 
show some correspondence with the discrete radio 
knots. Most soft X-ray emission, however, is from 
the disks and nuclear regions, rather than from the 
overlap region, even though several X-ray emission 
knots exist in the overlap. Yet, more soft X-ray 
structures in regions of the ISO "hotspots" are 



evidently better revealed in the new Chandra X- 
ray images. Both the X-ray and radio continuum 
show extended morphology (prominent ring struc- 
ture) in the northern galaxy NGC 4038, though 
not in NGC 4039. Since supernova remnants are 
most likely responsible for the majority of both ex- 
tended emission, these observations may be used 
as probes to identify past active star formation 
sites. According to the various models, the first 
burst of star formation may have occurred several 
~ 10 8 years ago when the first closest approach 
of the galaxy pair happened. Many X-ray knots 
with the radio emission correspondence can be 
even older, and they are reminiscent of the possi- 
ble star forming sites in the past, prior to or during 
the first strongest impact. Overall, these emission 
could mark roughly the sites of the interaction en- 
hanced large scale star formation across over the 
pair's disks at the earliest stage of merging. 

Ha emission and young star clusters shown in 
the HST images (Whitmore et al. 1999) indicate 
that recent vigorous star formation is also proceed- 
ing throughout the entire merging disks. Molec- 
ular gas concentrations in the overlap may have 
hidden some extremely young star clusters and HII 
regions. Nevertheless, these again are the rather 
recent large scale bursts of star formation occur- 
ring probably right after the genesis of the tidal 
tails, after the pair's first closest approach. 

Both the MIR and FIR observations (Mirabel 
et al. 1998; Xu et al. 2000; Evans et al. 1997; 
Bushouse et al. 1998) show the strongest peaks 
in the rather extended overlap region where the 
bulk of the molecular gas resides. But the MIR 
and FIR emission peaks differ from each other, 
whereas the FIR peaks roughly at same position 
as that of the CO. Generally speaking, the overlap 
region is the current sites of vigorous star forma- 
tion or starburst. Yet, the overall starburst level 
is at most modest. The strongest starburst sites 
with SFE comparable to that of LIGs/ULIGs are 
confined to ,6 kpc scale small regions in the over- 
lap, and the global SFE level is about same as 
those in the Milky Way disk GMCs. Despite some 
localized starbursts are still going-on, the recent 
peak of starburst phase in Arp 244 has most likely 
passed already, and the entire system is currently 
in a rather quiescent star formation phase. 

The radio continuum-to-CO ratio map implies 
low SFE at both nuclei (Fig. 7). The highest gas 
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concentration in Arp 244 revealed by the interfer- 
ometers is, however, in the northern nuclear re- 
gion, which has only moderate £ kpc scale cir- 
cumnuclear starburst ring, rather than a sub-kpc 
nuclear starburst. The highest radio-to-CO ratios 
are the confined bona fide starburst sites (Fig. 6) 
which appear to have some correspondence with 
the [C n] emission peaks (Nikola et al. 1998). 
The CO emission peaks in the overlap region ap- 
pear not to be at the exact locations of the local- 
ized starbursts, and the most intense starbursts 
might not be at an extremely high SFE level as 
in ULIGs. These highest molecular gas concentra- 
tions are likely to undergo future ultraluminous ex- 
treme starburst once they all merge together with 
the nuclear gas concentrations, and collapse into 
~ kpc scale in the final coalescence. Although the 
progenitor galaxies of Arp 244 could only have 
small stellar bulges, the bulges can be built up 
through merging (e.g., Carlberg 1999) in addition 
to nuclear gas infalls and concentrations. Thus, 
Arp 244 fits most likely between the two extremes 
modeled by Mihos & Hernquist (1996), and more 
than an order of magnitude enhanced SFE can be 
achieved in the final merging. 

Although most (> 60%) HI gas is distributed 
in the tidal tails (van der Hulst 1979; Hibbard 
ct al. 2000), there are still - 2 x 1O 9 M HI in 
the merging disks, whereas most ULIGs have lit- 
tle HI left in the merged disks (Hibbard & Yun 
1996, 2000; Mirabel & Sanders 1988, 1989). Most 
HI in Arp 244 is in the southern tail, where we 
tentatively detected CO emission at the tip of 
the tidal tail, indicating a molecular gas mass of 
^,2 x 10 8 M Q in the tidal dwarf galaxy. These are 
additional gas reservoir for future star formation, 
especially when the extraordinary HI tails rain 
back onto the merging disks (Hibbard et al. 1994), 
converting most HI gas into molecular phase. 

Our recent HI observations of NGC 6670 sug- 
gest that a precursor to Arp 244 can be that the HI 
disks have merged into a huge overlap concentra- 
tion in the extended HI disks, prior to the merging 
of the stellar and molecular gas disks (Wang et al. 
2000). In Arp 244, although an anti-correlation 
between the molecular and atomic gas distribu- 
tions appears to exist, there is still a significant 
amount of HI gas in the molecular gas overlap re- 
gion. Whether there was a HI gas overlap formed 
previously, prior to the merging of the molecu- 



lar gas disks, requires high resolution HI observa- 
tions and further comparative studies of the multi- 
wavelength observations, in combination with the 
numerical modelings, to test. Clearly, understand- 
ing of the formation and evolution of the overlap 
region is the key to ultimately comprehend the en- 
tire star formation history of Arp 244. 

4.4. Overlap Star Formation and Star- 
burst Mechanisms 

The formation of the overlap star formation re- 
gions, unnecessarily as overwhelmingly dominant 
as in the case of Arp 244, might be a quite common 
phenomenon in merging galaxies (Xu et al. 2000). 
Both the simulations and observations have shown 
that the gas is being transported into the nuclear 
regions when a pair of spiral galaxies undergoes 
through the merging process, in addition to be- 
ing dragged out into the tidal tails (e.g., Olson 
& Kwan 1990a,b; Noguchi 1991; Mihos & Hern- 
quist 1996; Barnes & Hernquist 1991, 1996; Scov- 
ille et al. 1991; Gao et al. 1999; Hibbard & van 
Gorkom 1996). Yet, models appear to be unable 
to reveal the early formation of the gas concen- 
tration in the overlap region when galaxies are 
still quite widely separated, like the intermedi- 
ate merger Arp 244 with the molecular gas over- 
lap, where most gas resides, and the early merger 
NGC 6670 with the atomic gas overlap. Could 
a significant amount of the atomic gas, especially 
in the southern progenitor, have ended up in the 
overlap region in the first place, prior to the forma- 
tion of the molecular gas overlap, as a result of ef- 
ficient HI cloud-cloud collisions than the collisions 
of the GMCs, when the two gas-rich spiral pro- 
genitors first collided? This is because the mean 
free path of a HI cloud is ~ 33 pc, much smaller 
than the size of the overlap region, yet the GMCs' 
mean free path is orders of magnitude larger and 
the chance of GMC collisions in the overlap is ex- 
tremely small (Jog & Solomon 1992). Just like 
the early merger NGC 6670, which has the HI gas 
concentration in the HI disk overlap region formed 
prior to the merging of the stellar and molecular 
gas disks (Wang et al. 2000), there might be a 
similar HI gas overlap formed in Arp 244 during 
its earliest merging stage. 

Jog & Solomon (1992) argued that a starburst 
occurs when the preexisting GMCs in the overlap 
region undergo radiative shock compressions by 
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the preexisting high pressure of the central molec- 
ular intercloud medium produced by the heating 
from the HI cloud-cloud collisions. If indeed an HI 
disk-disk overlap region could be formed prior to 
the merging of the stellar and molecular gas disks, 
the dominant star formation should not be occur- 
ring in the HI gas overlap region, given few GMCs 
exist in this HI cloud-cloud collision ISM during 
the early stage of merging. Instead, large-scale 
star formation and nuclear starbursts induced by 
the gas inflow are probably at play (Combes et al. 
1994). It is later formation of the molecular gas 
disk overlap, within this preexisting high pressure 
overlapping HI gas concentration, when the merg- 
ing advances into the intermediate stage, that will 
probably make GMCs in the overlap undergo star- 
burst by the radiative shock compressions (Jog & 
Solomon 1992). 

Although the soft X-ray emission in the overlap 
region in Arp 244 is neither strong nor extended, 
as revealed by the ROSAT HRI (Fabbiano et al. 
1997), this may be due to the high gas column 
density in the overlap which absorbs the most ex- 
tended soft X-ray emission. The Chandra soft X- 
ray images of much improved sensitivity and reso- 
lution do show some extended features in the over- 
lap, which are particularly in good correspondence 
with the MIR "hotspots". This is in the right 
track, just as expected from the prediction of the 
overlap starburst model of Jog & Solomon (1992). 
Higher energy X-ray imaging that can penetrate 
through the heavy dust/gas concentration in the 
overlap, yet still orders of magnitude below the 
Compton-thick limit, and the extinction-corrected 
soft X-ray images of Arp 244 could be the key 
to further test the overlap overpressure starburst 
mechanism. 

In advance mergers, the double-nucleus gas 
disks are already merging with the overlap gas con- 
centrations, with most gas either between (e.g., 
NGC 6240, Tacconi et al. 1999) or around the 
double-nucleus (e.g., Arp 220, Downes & Solomon; 
Sakamoto et al. 1999) in ^kpc scale, in addition 
to the nuclear gas concentrations. Thus, a possible 
merging sequence for a gas-rich pair of spirals may 
look like the following: from the build-up of the HI 
gas disk overlap regions (NGC 6670) — > the suf- 
ficient overlap in both stellar and molecular disks 
(Arp 244) — > the kpc scale merging double-nucleus 
gas disks, either with a central gas concentration 



between the nuclei (NGC 6240), or with an ex- 
tended gas disk of several kpc surrounding the nu- 
clei (Arp 220). The transformation of the merging 
phases and the further development of the over- 
lap regions are also accompanied by the change 
in the dominant sources of the total power out- 
put: from large scale star formation within galaxy 
disks in early mergers — > the mainly localized star- 
bursting overlap regions between the galaxies in 
intermediate mergers — > the extreme starburst in 
the highly concentrated double-nucleus sub-kpc 
sources in advanced ULIGs. 

The overpressure starburst mechanism (Jog & 
Solomon 1992) is probably happening in interme- 
diate mergers, especially those with the extensive 
overlap regions like Arp 244. Yet, this mechanism 
appears to require a sufficient overlap between the 
gas disks so that the early formation of the gas 
overlap region is possible. Is this the sole star- 
burst mechanism during the entire merging pro- 
cess when two gas-rich spiral galaxies merge? It 
is likely that different starburst mechanisms play 
various roles during different merging stages along 
a merger sequence. What then are the other trig- 
gering mechanisms for the starbursts? Wilson 
et al. (2000) compared their new OVRO CO map 
with the ISO MIR emission and suggested that 
molecular cloud collisions may play an important 
role for the local intense starburst and the strong 
MIR emission. Yet, Lo et al. (2000) found no 
obvious signatures of cloud collisions from ana- 
lyzing the molecular gas kinematics in the BIMA 
+ the NRAO 12m full-synthesis CO data. Per- 
haps Arp 244 is on its way of changing from the 
overpressure overlap starburst to the direct cloud- 
cloud collision starburst. 

An extensive active star forming molecular gas 
overlap region observed in Arp 244 may not al- 
ways be the case for most mergers. In fact, most 
LIGs do not show such an extensive molecular gas 
overlap region (e.g., Gao et al. 1997, 1999; Bryant 
& Scoville 1999). Even though star formation in 
the overlap can be common in mergers, the over- 
lap starburst is probably much less dramatic than 
that of Arp 244 in most cases (e.g., Xu et al. 2000). 
But, do all gas- rich mergers develop such an exten- 
sive molecular gas overlap phase during the course 
of merging? We do not have a definite answer. 
Only in late merging stage have both observations 
and numerical simulations demonstrated the exis- 
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tence of the highly concentrated molecular gas be- 
tween the merging nuclei. Also, only in advanced 
mergers, can the chance of GMC-GMC collisions 
be much enhanced since almost all the molecu- 
lar gas (in both the overlap and the two nuclear 
gas concentrations) finally condenses into a kilo- 
parsec scale. Now, the mean free path of GMCs 
is less than the size scale of the gas concentra- 
tion region since the volume filling factor of GMCs 
is ten-fold increased, making even the intercloud 
medium molecular. Therefore, direct collisions of 
GMCs could be playing an important role for the 
enhanced starburst when merging progresses into 
advanced stage in the final coalescence. 

5. CONCLUSIONS 

We summarize our main results and present our 
concluding remarks in the following: 

1. Our fully sampled CO (1-0) map (FWHM~ 
55") at half-beam spacing reveals a factor of 5.7 
larger total CO flux in Arp 244 than the previously 
accepted. This is mainly due to the enormous CO 
extent of Arp 244, which is much larger than the 
single pointing 55" beam used in the old obser- 
vation. In general, single beam observations of 
merging systems may well systematically underes- 
timate CO content in the cases where the merg- 
ing systems are much larger than the single beam, 
unless a map of extensive coverage of the merging 
systems is made. Most of CO emission in Arp 244 
indeed peaks at the disk overlapping region with 
a broad velocity spread of 600 kms -1 . We ob- 
tain a total molecular gas mass of 1.5xl0 10 M Q 
for Arp 244 using a standard CO-to-H 2 conversion 
factor. 

2. HCN(l-O) emission has been detected in 
Arp 244. HCN observations suggest that there 
is only a small amount of dense molecular gas in 
Arp 244 compared to the total molecular gas con- 
tent. The fraction of the dense molecular gas is 
only comparable to that of normal spiral galax- 
ies and is much lower than that of luminous and 
ultraluminous infrared galaxies (LIGs/ULIGs). 

3. We detected extended CO emission far away 
(nearly two full beams) from the major CO con- 
centrations, such as those revealed by the inter- 
ferometers. The peak CO emission in the over- 
lap region has an average gas surface density near 
~ 200 M Q pc -2 over ~ 5 kpc — the scale probed 



by the 55" beam. Even the resolved interferometry 
CO images only reveal that the highest gas surface 
density is about 10 3 M Q pc~ 2 , still order of magni- 
tude lower than that of ULIGs. This is consistent 
with the small fraction of dense gas found from 
the HCN measurements. Weak CO emission from 
the tidal dwarf galaxy at the tip of the southern 
tail has also been tentatively detected and the es- 
timated molecular gas mass is possibly larger than 
- 2 x 10 8 M Q . 

4. The excellent correspondence of the radio 
continuum images with the FIR maps is evident 
in Arp 244, which also implies the tight correla- 
tion between the two emission. We use the high 
resolution VLA 20cm image to approximately rep- 
resent the FIR emission distribution in Arp 244, 
and compare with the CO images of comparable 
resolution to estimate the local star formation ef- 
ficiency (SFE), which is defined as the FIR to 
molecular gas mass ratio, and is now indicated 
by the radio-to-CO ratio. SFE is low throughout 
the entire system with a global SFE= 4.2 Lq/Mq, 
just same as that of the Galactic disk GMCs. Only 
some confined regions in the overlap, inconspicu- 
ous in optical observations, are the sites of the cur- 
rent powerful starbursts with SFE >,20 L Q /M©, 
comparable to that of LIGs/ULIGs. The VLA ra- 
dio continuum images also reveal a very similar 
morphology as that of the molecular gas distribu- 
tion. 

5. Although some localized starbursts are going 
on in the confined regions in the overlap, the bulk 
of the molecular gas is forming into stars with a 
low SFE, just same as that of normal spiral galax- 
ies. The large amount of total molecular gas, yet 
only a small amount of dense molecular gas and 
a low SFE, revealed from our observations, com- 
bined with the theoretical predictions of models 
of gas-rich mergers, seem to suggest that Arp 244 
has the ultraluminous extreme starburst potential 
when merging proceeds into late stage. Arp 244 
appears to be a snapshot of an ULIG in its early 
stage of making. 

6. Multi- wavelength comparison together with 
numerical models appear to indicate that Arp 244 
has recently passed the peak of large scale active 
star formation. Globally, the bulk of the molec- 
ular gas in Arp 244 is currently not in a star- 
burst phase. The confined bona fide starburst 
sites of the highest SFE appear to be offsets from 
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the peaks of the CO, radio continuum, FIR, and 
MIR emission. The future ultraluminous extreme 
starburst phase can possibly be reached once the 
CO concentrations in the overlap merge with the 
double-nucleus gas concentrations in the final co- 
alescence. 

7. In conclusion, our newly determined 5.7 
times larger total molecular gas mass in the "An- 
tennae" galaxies lowers the global SFE to be com- 
parable with that of GMCs in the Milky Way 
disk. Comparing the CO images with other ob- 
servations, including our own VLA 20 cm con- 
tinuum imaging, an overall relatively low molec- 
ular gas density and a low SFE across over the 
bulk of the molecular gas distribution have been 
revealed, which is in sharp contrast with most 
LIGs/ULIGs. The strongest starburst sites are 
currently confined in small regions heavily ob- 
scured in the dusty patches. We conclude that 
the large reservoir of molecular gas is the sufficient 
fuel for Arp 244 to ultimately enter an ultralumi- 
nous extreme starburst phase as merging proceeds 
to the advanced stage. The ultraluminous phase 
is possibly achieved once almost all molecular gas 
condenses into ^kpc scale, with an order of magni- 
tude increase in the gas surface density and more 
than ten-fold higher SFE. 
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Table 1 

Summary of the NRAO 12m CO Mapping Results 



No. R.A. a 


Decl. a 


/co b 


V co 


AVfwhm 


arcsec 


arcsec 


K kms- 1 


kms -1 


kms -1 



1 


0.0 


0.0 


19.5 


1632 


93 


2 


0.0 


27.5 


11.9 


1612 


204 


3 


0.0 


-27.5 


20.1 


1607 


201 


4 


0.0 


55.0 


5.2 


1607 


282 


5 


0.0 


82.5 


3.5 


1561 


331 


6 


0.0 


-82.5 


10.6 


1659 


281 


7 


0.0 


-55.0 


19.2 


1604 


227 


8 


0.0 


110.0 


0.4 C 


1663 


110 


9 


0.0 


-110.0 


3.9 


1638 


348 


10 


0.0 


-137.5 


1.3 


1660 


229 


11 


27.0 


0.0 


21.5 


1593 


198 


12 


27.0 


27.5 


13.6 


1610 


280 


13 


27.0 


-27.5 


26.9 


1565 


200 


14 


27.0 


55.0 


3.4 


1567 


229 


15 


27.0 


-55.0 


27.4 


1567 


245 


16 


27.0 


82.5 


2.5 


1553 


289 


17 


27.0 


-82.5 


14.0 


1604 


270 


18 


27.0 


110.0 


d 






19 


27.0 


-110.0 


2.6 


1698 


166 


20 


27.0 


-137.5 


1.2 C 


1618 


308 


21 


-27.0 


0.0 


14.4 


1627 


165 


22 


-27.0 


27.5 


9.6 


1609 


144 


23 


-27.0 


-27.5 


13.0 


1598 


255 


24 


-27.0 


55.0 


4.4 


1582 


150 


25 


-27.0 


-55.0 


7.7 


1621 


279 


26 


-27.0 


82.5 


2.7 


1626 


241 


27 


-27.0 


-82.5 


5.3 


1658 


321 


28 


-27.0 


110.0 


2.2 


1568 


463 


29 


-27.0 


-110.0 


3.7 


1640 


302 


30 


-27.0 


-137.5 


d 






31 


53.9 


0.0 


11.8 


1592 


213 


32 


53.9 


27.5 


7.7 


1572 


246 


33 


53.9 


-27.5 


18.9 


1558 


215 


34 


53.9 


55.0 


2.4 


1648 


204 


35 


53.9 


-55.0 


19.8 


1560 


193 


36 


53.9 


82.5 


1.2 


1652 


241 


37 


53.9 


-82.5 


9.1 


1597 


254 


38 


53.9 


-110.0 


2.7 


1630 


360 


39 


53.9 


-137.5 


1.5 


1651 


283 


40 


-53.9 


0.0 


9.4 


1628 


241 


41 


-53.9 


27.5 


3.3 


1622 


222 
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Table 1 — Continued 



No. 


R.A. a 


Dccl. a 


/co b 


^co 


AVfwhm 




arcsec 


arcsec 


K kms- 1 


kms -1 


kms -1 


A 1 

42 


con 
— 53.9 


ri7 c 

— 27.5 


A A 

4.4 


1669 


1 A C 

245 


A O 

43 


con 
— 53.9 


c c n 
55. L) 


1 1 

2.2 


1616 


o o o 
388 


A A 

44 


con 
— 53.9 


c c n 
— 55.0 


1 O 

2.3 


1696 


£i 
02 


/I c 

45 


con 
— 53.9 


o o c 
82.5 


1 1 c 

1.1 


1 C /I 7 

I54r 


1 c n 
15U 


46 


con 
— 53.9 


on c 

— 82.5 


n nc 
0.9 


1700 

1 (88 


1 c 

25 


A 7 

4 f 


co n 
— 5o.9 


1 1 n n 
— 11U.U 


1 1 

z.z 


i cno 
1595 


OIO 

SIS 


45 


on n 

81). 9 


n n 
U.U 


/I c 
4.0 


10U ( 


O 1 £ 

olO 


49 


on n 
51). 9 


07 c 

z 1 .5 


1 
Z.J 


1 C£7 

loO I 


IOC 

zo5 


5U 


on n 
51). 9 


07 c 

— z / .5 


R C 

0.5 


1 r/J7 

150 ( 


110 

Zz5 


C 1 

51 


on n 
51). 9 


c c n 
55. U 


n £c 
U.O 


1 70 /I 

1 ( o4 


1 A 7 

14/ 


CO 
5z 


on n 
51). 9 


1 07 C 

— loV .5 


1 1C 

l.z 


1 71 C 

1 ( z5 


1 17 
IZY 


CO 

5o 


on n 
— 5U.9 


n n 
U.U 


1 7 
1. / 


104Z 


CO 

oS 


c /i 
54 


on n 
— 5U.9 


07 C 

ZY.5 


1 /I 

Z.4 


105U 


ono 
SUS 


r r 

55 


on n 
— 8U.9 


17 C 

— 27.5 


1.0 


1632 


i /) i 
242 


c 

56 


on n 
— 8U.9 


c c n 
55. U 


n ac 
U.O 


1 £;nn 
loUU 


1 C 7 

157 


C 7 

57 


on n 
— 8U.9 


c c n 
— 55. U 


n /i c 
U.4 


1 a o o 
1638 


1 1 
21 


c o 
5o 


on n 
— 8U.9 


on c 

— 82.5 


n cd 
U.5 


1 C 71 

1571 


7 /I 

<4 


c n 
59 


i nn o 
1U9.3 


n n 
U.U 


i i 
1.2 


1554 


18b 


£n 

60 


i nn o 
1U9.3 


17 C 

27.5 


1 oc 

1.3 


i £in7 
169 I 


1 C 1 

152 




i nn o 
1U9.3 


17 C 

— 2/. 5 


n nc 
U.9 


lolU 


1 7n 

179 


a i 
62 


1 (17 n 
— 1U7.9 


n n 
U.U 


d 






a o 
63 


i (17 n 
— 1U7.9 


17 C 

27.5 


1 nc 
l.U 


i c ni 
15U2 


1 nn 
1UU 


64 


-107.9 


-27.5 


0.4 C 


1599 


54 


65° 


81.8 


-55.0 


4.6 


1538 


213 


66° 


109.3 


-82.5 


1.7 


1584 


178 


67° 


81.8 


-82.5 


3.5 


1575 


284 


68° 


81.8 


-110.0 


0.9 


1629 


251 


69° 


109.3 


-55.0 


1.7 


1564 


201 


70° 


109.3 


-110.0 


1.7 


1617 


324 


71° 


136.8 


-55.0 


0.4 C 


1588 


30 


72° 


136.8 


-82.5 


1.2 C 


1644 


370 


73° 


136.8 


-110.0 


0.5 d 


1520 


300 


74 


25.9 


-39.1 


34.3 


1558 


209 



a Offsets relative to RA=12' l 01 m 53.1 s , Dcc=-18°52'05."0 (J2000, the nucleus of NGC 4038). 

integrated CO line intensity, which is calculated using the same broad line emission window of 1300 to 
1900 kms -1 . The Jy/K (T£) conversion adopted is ~ 35 Jy/K. Some offsets have slightly larger Ico when 
a narrower line emission window is used. 

c Tentative detections (^3cr). 

d Non-detections (< 3er). 

e Deep CO spectra are shown in Figure 5a. 
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Table 2 

Global Properties of Arp 244 



Parameter 


Value 


References 


dh 


20 Mpc 


van dcr Hulst 1979; Mirabel et al. 1998 


Lb 


2.9 x 1O 1O L 


from B mag. in RC3 


Lm 


6.2 x 1O 1O L (IRAS) 


Soifer et al. 1989; Sanders ct al. 2001 in prep. 


Lm 


8.0 x 1O 1O L (KAO) 
90.6 K km s- 1 (T*) 
3172Jykms- 1 


Bushouse et al. 1998 


Ico 


this work 


ScodV 


this work 


Leo 


0.3 x 10 10 K kms- 1 pc 2 (T mb ) 


this work 


M(H 2 ) 


1.5xl0 10 M Q 


this work 


£ m /M(H 2 ) 


4.2 L /M 


this work 


M(HI) 


0.4 x 1O 1O M 


van der Hulst 1979; Hibbard et al. 2001 in prep. 


M dus t 


- 10 8 M 


Haas ct al. 2000 


M gas /M dust 


- 190 


this work 


Lhcn 


0.7 x 10 8 K kms" 1 pc 2 (T mb ) 


this work 


Lhcn/ Leo 


0.02 


this work 


M(H 2 )tidal dwarf 


>2 x 10 8 M 


this work 



0010128.figl.gif 

Fig. 1. DSS image of Arp 244 (the "Antennae" galaxies) and the locations of our deep HCN and 
CO search. Thin circles are CO beams, observed at southeastern edge of the disk overlap where the southern 
tidal tail begins, and at the tidal dwarf galaxy near the tip of the southern tail. The thick circles mark the 
two HCN beams. The southern HCN beam covers both the overlap and the nuclear region of NGC 4039. 
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0010128.fig2.gif 

Fig. 2. Fully mapped NRAO 12m CO (1-0) integrated intensity map (in contours) overlaid on the 
HST/WFPC2 image. The contours are 3, 4, 6 to 28 K kms -1 (in increases of 2 K kins" 1 on a scale). 
The CO spectra at each of the observed position are also plotted (velocity ranges from 1150 to 2050 kms -1 , 
and TjJ ranges from —25 to 155 mK). The currently accepted total molecular gas mass in the literature was 
based upon only one beam measurement near the CO peak position (Sanders & Mirabel 1985). Obviously 
one NRAO 12m beam (top left) does not capture anywhere even near the CO emission from the overlap 
region. We have sampled 73 positions at half-beam spacing to ensure that no further extended CO emission 
was significantly detected in the outer edge, in typically 1.5 hr of the integration time. 
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Fig. 3. 

CO spectra closest to the peak CO emission in the overlap region compared with the old CO spectrum of 
Sanders & Mirabel (1985). We detected broader CO profile, which was previously missed due to insufficient 
velocity coverage in the old CO observation, and stronger CO emission. 
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0010128.fig4a-4e.gif 

Fig. 4a — 4e. CO integrated intensity maps (contours) in 100 km s^ 1 velocity width. All contours plotted are 
0.4, 0.7 to 2.5 K kms -1 in increases of 0.3 K kms" 1 . 
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0010128.fig4f.gif 

Fig. 4f. The CO integrated intensity map over the total velocity range (essentially same as the con- 
tour map in Figure 2, but here all CO data have been used), which has contours 1.4, 2, 3, 4, 6 to 34 
K kins" 1 , in increases of 2 K kms" 1 . 
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Fig. 5a. The 

deep CO(l-O) spectra near the southeastern (SE) edge of the overlap region (the SE corner of the merging 
disks in Fig. 1). The numbers refer to the offsets in arcsec, with respect to the northern nucleus, and the 

central position here corresponds to an offset of (110", —82.5"). Fig. 5b. Both the summed average 
spectrum (top) from all 5 positions searched for CO (more than 27 hours integration) in the tidal dwarf 
galaxy, and the deep integration spectrum (bottom) at position of RA=12:0f :26.0, DEC=— 19:00:37.5 
(J2000) (more than 12 hour integration), show tentative weak detections. Fig. 5c. HCN(l-O) spectra 
obtained from the northern galaxy NGC 4038 anc^from the overlap region plus the southern galaxy 

NGC 4039. 
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Fig. 5c 
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0010128.fig6.gif 

Fig. 6. The 20cm radio continuum to CO(l-O) emission ratio map (contours) compared with the 
HST/WFPC2 image (gray-scale). This illustrates the sites of the most intense starbursts of the highest 
star formation efficiency (SFE). White contours indicate the sites of the highest ratios, ~ 5 times above the 
average, suggesting the highest SFE ^,20 L Q /M Q . Black contours are the sites of a factor of 2 lower in the 
radio-to-CO ratios, implying a SFE ^,10 L Q /M0. 
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0010128.fig7.gif 

Fig. 7. The overall ratio map (contours) of the 20cm radio continuum to the CO emission compared 
with the VLA 20cm radio continuum image (gray-scale). Here black contours are the sites of the highest 
radio-to-CO ratios (same as those white contours in Figure 6) and the ones a factor of 1.5 lower, gray 
contours are for the radio-to-CO ratios of a factor of 2 lower than the highest (same as those black contours 
in Figure 6) and those a factor of 1.5 further lower, and white contours here are for the lowest radio-to-CO 
ratios, a factor of 4 and 6 lower than the highest, indicating the sites of SFE ~ 5L /M and lower, about 
the average SFE across over the entire system. 
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